Absence of an isotope effect in the magnetic resonance in high-T c superconductors 
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An inelastic neutron scattering experiment has been performed in the high-temperature supercon- 
ductor YBa2Cu3 06.89 to search for an oxygen-isotope shift of the well-known magnetic resonance 
mode at 41 meV. Contrary to a recent prediction (I. Eremin, et al., Phys. Rev. B 69, 094517 
(2004)), a negligible shift (at best < +0.2 meV) of the resonance energy is observed upon oxygen 
isotope substitution ( 16 0— > ls O). This suggests a negligible spin-phonon interaction in the high-T c 
cuprates at optimal doping. 



In conventional superconductors, pairing between elec- 
trons is mediated by lattice vibrational. This has been 
demonstrated by an isotope effect on the superconduct- 
ing (SC) transition temperature, T c . In high-Tc copper 
oxides superconductors T c exhibits a weak shift at op- 
timal doping upon isotope substitution 2 which increases 
at lower doping. In particular, the oxygen-isotope shift 
( 16 0^ 18 0) has been extensively studied 3-4-5 . At opti- 
mal doping, a small isotope-effect exponent is deduced 
otT c = —d\nT c /dh\M ~ 0.05 much lower than the | 
value expected from pure electron-phonon interaction, 
casting some doubt on a superconducting mechanism me- 
diated by phonons. Further, the proximity of the antifer- 
romagnetic (AF) insulating state and the unconventional 
d-wave symmetry of the SC gap favored mechanisms for 
high T c superconductivity where electron-electron (el- 
el) interactions predominate. However, there has been 
a revival of interest in electron-phonon coupling as sev- 
eral experiments point towards a non-negligible electron- 
phonon interaction^ 7 ^. In particular, the "kink" change 
of slope of electronic dispersion observed around ~ 
70 meV by angle-resolved photoemission spectroscopy 
(ARPES) along the nodal direction in various cuprates 
could be interpreted as an electronic coupling to a phonon 
mode 6 . As a matter of fact, various physical properties 
such as penetration depth? or ARPES spectral display 
relatively large isotope effects thus highlighting the open 
question: what can be the role of phonons in determining 
the superconducting properties of cuprates? 

On the other hand, the spin excitation spectrum of 
the copper oxide superconductors is particularly rich. 
Above T c , magnetic fluctuations are mainly observed 
around the planar wavevector Qaf = (7r/a, 7r/a) char- 
acteristic of antiferromagnetism (AF) in the undoped 
parent compoundsSiiSiii*^. Below T c , a collective mag- 
netic mode, referred to as the "resonance peak", ap- 
pears at a well-defined energjiSiiSiiiiiSii 3 ^ at (ir/a, tt/o) 
and exhibits strong dispersions for wavevectors around 
Qa i 15 i 16 i 17 . This mode is now observed in all high- 
Tc superconductors systems studied by inelastic neutron 
scattering (INS) experiments2iiS*iS whose maximum T c 



reaches 90 K. Recently, an analogous feature has even 
been reported in the single-layer material La2- x Sr x Cu04 
as wellSS. Depending, or not, whether the magnetic fluc- 
tuations are observed in the normal state, the resonance 
peak either corresponds to a modification of magnetic 
spectrum in both momentum and energy (in underdoped 



cuprates) 



or simply emerges from the magnetic 



electron-hole continuum (in optimally doped and over- 
doped cuprates) 14 i 17 i 21 . 

This mode is typically assigned to an excitonic 
bound state in the superconductivity-induced gap 
in the spectrum of electron-hole spin-flip Stoner 
excitationa 22 ! 23 ! 24 ! 25 ! 26 ! 27 ! 28 ! 29 ! 30 . Within that frame- 
work, the resonant mode is a direct consequence of un- 
conventional superconductivity of cJ-wave symmetry oc- 
curring in the high-T c copper oxides. This approach 
is particularly suited for optimally doped superconduc- 
tors where normal-state magnetic fluctuations are con- 
sistent with a broad magnetic electron-hole continuum. 
Further, the most recent neutron developments 1 ^ allow 
"silent bands" (where the magnetic collective mode is 
overdamped) to be related to the detailed momentum 
shape of the Stoner continuum expected from the Fermi 
surface topology and the d-wave superconducting order 
parameter. This connection has been explicitly assigned 
in ref^i. This spin exciton mode can be derived from 
an effective t — J HamiltonianSSiSS. Within that model, 
an interesting proposal has been made^S that magnetic 
properties could display significant isotopic effects if both 
the hopping integral, t, and the superexchange interac- 
tion between neighboring spins, J, are renormalized by 
phonons. By changing the oxygen isotope 16 by 18 0, 
they predict a change of the resonance peak position of a 
few meV, mostly due to a re-normalization of the hopping 
integral. 

Thus, by measuring the isotope dependence of the po- 
sition of the resonance peak, INS could provide direct ev- 
idence for the presence of the electron-phonon coupling 
in cuprates. In contrast, we here report the probable ab- 
sence of an isotope effect in the magnetic resonance peak 
in YBa 2 Cu 3 6 . 8 9 (YBCO). 
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FIG. 1: (Color online) Magnetic susceptibility of both sam- 
ple arrays used in this experiment. Onset superconducting 
temperatures occur at 88.8 K for samples with 18 and 90 
K for 16 0. The curve of each mounting has been obtained 
by weighting the susceptibility curve of each individual single 
crystal by their mass. The dispersion of T c for the samples of 
each array is less than 0.3 K. 




35 40 45 
Energy (meV) 



55 



FIG. 2: (Color online) Neutron intensity measured at Q=(- 
0.5,-0.5,10.3) and at T=100 K showing an oxygen phonon 
mode in both 16 O sample and in the 18 O sample. 



The inelastic neutron experiment was performed on the 
IT spectrometer at Laboratoire Leon Brillouin (LLB) in 
Saclay. The spectrometer used a vertically and horizon- 
tally focusing monochromator and analyzer, comprising 
Cu (111) crystals and pyrolytic graphite (PG002) crys- 
tals, respectively. The measurements were performed 
with a fixed final neutron energy of 30.5 meV. A filter 
was inserted into the scattered beam in order to elim- 
inate higher order contamination. The crystals were 
oriented such that momentum transfers Q of the form 
Q = (H, H, L) were accessible. We use a notation in 
which Q is indexed in units of the tetragonal reciprocal 
lattice vectors 2n/a = 1.63A -1 and 2-njc = 0.54A -1 . 

High-quality single crystals of YBa2Cu3 6 07_5, of typ- 
ical mass ~ 0.1-0.2 g, were prepared. They have been 
separated in two distinct sets of similar total mass (~ 
0.6 g). The two batches of samples were mounted in a 
furnace in separate quartz tubes side by side to ensure 
identical thermal history. One tube was charged with 
high purity 16 O oxygen gas while the other was charged 
with 99% 18 O enriched oxygen. After annealing for 24 
hours at 830 °C the pair of tubes was removed, evacuated, 
recharged and reannealed for a total of 10 exchanges to 
ensure maximal isotope exchange. Both sample sets were 
then slow cooled over 48 hours to 550 °C then annealed 
there for 10 days to ensure the same uniform oxygena- 
tion of the chains. The samples were lightly underdoped 
as confirmed by thermoelectric power measurements 31 . 
Based on previous annealing experience for YBCO ce- 
ramics and crystals^ we expect the samples to have an 
oxygen content of w 6.89 and an underdoped hole con- 
centration of rih ~ 0.145. The magnetic susceptibility for 
each individual sample has been measured by a SQUID 
( superconducting quantum interference device) magne- 



tometer. The crystals for each individual sample were 
then aligned on two distinct arrays (referred hereafter as 
Si6 and Sis) of similar volume, each array being made 
of about four single crystals. The magnetic susceptibil- 
ity measurements of each sample mounting are shown 
on Fig. 2] one can identify the onset of the supercon- 
ducting transition at 90 K for Sie and 88.8 K for Sis- 
A difference in T c of ~ -1.2 K is then observed between 
both samples. It is actually larger than the reduction ex- 
pected from the usual isotope effect which is about -0.77 
± 0.2 K for similarly-doped YBCO 4 . The slight excess 
in isotopic shift in T c is possibly significant and may sug- 
gest a slightly lower oxygen content and doping state for 
the 18 O sample, despite the attempt to ensure identi- 
cal thermal histories. To address this issue, we further 
determined the c-lattice parameter of both samples. Ac- 
curate measurements using the triple axis 4F1 spectrom- 
eter installed on a cold source at LLB yield c= 11.674 A± 
0.004 for Si 6 and c= 11.678 A± 0.004 for S 18 . These c- 
lattice parameters are consistent with an oxygen content 
of about x=0.8922iS, and given the slope dc/dx = —0.11 
A22ii£, indicates that the 18 O sample has a possible lower 
oxygen content by Ax = 0.02 ±0.03. Thus, it is in agree- 
ment with the possibility that part of the difference in T c 
for each sample is related to a slight difference in doping. 
Using the relationship between oxygen content and the 
hole doping^ 4 ., this would correspond to a difference of ~ 
671^=0.004 ± 0.006 in doping. 

In order to check the isotope exchange process, we per- 
formed Raman scattering at room temperature as well as 
INS measurements of a particular oxygen phonon mode. 
The three oxygen modes observed in Raman scattering, 
namely the c-axis vibration of the apical oxygen and the 
in-phase and out-of-phase oxygen vibrations in the Cu02 
plane, revealed isotopic shifts to lower energy in 18 sam- 
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FIG. 3: (Color online) Resonant magnetic peak intensity in 
both 16 samples and 18 O samples: a) Raw neutron intensity 
at 12 K and 100 K for both samples, b) Difference of the 
neutron intensity 12 K - 100 K for each sample. The total 
counting time reached 1.5 hour per point at each temperature 
to reduce the error bar on the energy of the resonance to 0.12 
meV. The curves are not normalized by samples mass. 



pie with respect to the 16 sample. Assuming that they 
are pure oxygen modes, the fraction of exchanged oxy- 
gen is y~ 0.95 ± 0.05. Being a surface-sensitive tech- 
nique, Raman scattering does not indicate if the isotope 
exchange occurred within the bulk of the samples. This 
might be problematic as oxygen diffusion is very slow. 
We then measured with inelastic neutron scattering a 
particular phonon mode whose eigenvector is predomi- 
nantly related to a vibration along the c axis of the oxy- 
gens of the Cu02 planed. This phonon mode, measured 
at Q=(-0.5,-0.5,10.3) where its structure factor is larger, 
is found at 43.6 ± 0.1 meV in S i6 and 41.3 ± 0.1 meV in 
Sis (Fig. |2J). Again assuming that the phonon is a pure 
oxygen mode, this corresponds to an isotope exchange 
of r/=0.92 ± 0.08. In fact, the eigenvector for that spe- 
cific mode corresponds to about 90 % weighting by the 
oxygen atoms. Therefore, the energy shift of the phonon 
measured in INS is fully consistent with the nearly full 
oxygen exchange deduced from the Raman data. The 
oxygen exchange thus occurred throughout the bulk of 
the material thus confirming the sample homogeneity. 

We now turn to the central result of this report. Fol- 
lowing previous studiesi2*iiiSiii, we identify the reso- 
nant magnetic mode by constructing, for each sample, 
the difference between constant- Q scans measured at 12 
K (< T c ) and 100 K (> T c ) and at the wave vector Q=(- 
0.5,-0.5,5.1). The magnetic resonant mode in both sam- 
ples is shown in Fig. |3Jb as well as the raw energy scans 
for both samples in Fig. |3Ja. At T=100 K, the back- 



ground, which displays the same phonon mode as the 
mode shown in Fig. [21 but with much weaker intensity 13 , 
is also shifted upon the isotope substitution. This phonon 
is known to exhibit no temperature dependence across 
Tc- Further, the phonon scattering structure factor for 
the chosen wavevector with L=h. \ is reduced enough that 
the background subtraction procedure in the determina- 
tion of the resonance peak energy, its broadening and 
its amplitude does not affect significantly the observed 
difference of Fig. [2Jb. The resonance peak energy is 
found at u R = 41.26± 0.12 meV in S w and lo r = 41.31± 
0.12 meV in Sis- The two energies are therefore not 
distinguishable within errors. The two peaks exhibit a 
slight difference in amplitude as well as in width (7.1 ± 
0.3 meV for Si6 and and 7.8 ± 0.3 meV for Sis). In 
principle, the observed width is not intrinsic but is con- 
trolled by the convolution product of dispersive excita- 
tions around (tt, tt) with the spectrometer resolutioni^*^. 
The product of the peak amplitude by its width in Fig. 
Elb, representing the magnetic resonant spectral weight 
at (tt, 7r), is similar in both samples within errors. Us- 
ing the spectral weight of the phonon presented in Fig. 
13 one can calibrate the absolute magnetic intensity of 
the resonance pea k 11 ' 13 . For both samples, we deduce an 
energy-integrated magnetic spectral weight of 2.6±0.4 fi 2 B 
at the (tt, tt) wavevector, or 0.06±0.01 fi B for energy- and 
q-intcgrated magnetic spectral weight in agreement with 
a previous report for a similar doping levelii. 

We then basically observe no isotope effect of the mag- 
netic resonance peak: a ulR = — dlnun/dlnM ~ 0. To be 
complete, there is however the possibility of a slight dif- 
ference in doping between the two samples which might 
induce a slight change in the resonance energy if the res- 
onance energy is proportional to T c , as it is typically 
observecUfl*ii. (To what accuracy this proportionality 
strictly applies is still an open question). According to 
this empiric relation, the resonance peak energy could be 
renormalized by about 0.5% in the 18 sample as com- 
pared with the 16 sample, i.e. an energy shift of ~ 
-0.18 meV. In such a case, an isotope effect can be esti- 
mated of Slur ~ +0.23, yielding an isotope exponent of 
&res = —0.05. Therefore, the isotope shift of the reso- 
nance peak energy can be at most Sur < 0.23 meV ± 0.2 
meV. The deduced isotope-effect on the resonance peak is 
then very small and actually similar in magnitude to the 
small isotope effect of the superconducting transition for 
optimally doped cuprates, aj c = 0.05 2 , although with an 
opposite sign. It should be noticed that the overall effect 
might be simply overshoot by the difference in doping 
between both samples. 

The absence of an isotope effect on the resonance en- 
ergy is actually quite surprising as, in the spin exciton 
model, the bound state energy is very sensitive to both 
band structure, via the hopping integrals t,t' and the 
interactions, g. Within a random-phase approximation 
(RPA) scheme, the resonance energy at the AF wavevec- 
tor is usually defined as the pole of the interacting suscep- 
tibihty, 1 - g/2Re Xo (QAF ,ur) = O 23 . 24 ! 25 . 2 ^^ 2 ^ , 
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where Xo(Q, w ) K 1/* ^ s * ne ^ are s P m susceptibility of 
a (i-wave superconductor. The interaction g can be ei- 
ther the superexchange interactions^i^iSSiSSiSS, 4 Jaf at 
(tt, 7r), or some spin-fermion coupling 25 27-30 . Looking in 
more detail at the pole condition, one clearly sees that a 
shift in g and in t have opposite effects on the resonance 
energy ujji: if the interaction is reduced, the resonance 
energy will increase whereas if the band hopping inte- 
gral is reduced the resonance energy will decrease. Using 
such a model, Eremin et al2£ expected a shift of about -2 
meV of the resonance peak from 16 O to 18 O, correspond- 
ing to a UR — 0.4. That was because the band structure 
hopping integral t is thought to be the quantity strongly 
dependent on the electron-phonon coupling. As the ob- 
served 8 UR has an opposite sign (if any) of the theoreti- 
cal prediction 2 ^, it is doubtful that the electron-phonon 
coupling is renormalizing the band structure hopping in- 
tegral. Within that model, one can nevertheless explain 
the observed positive sign of 5 UJR by a small renormal- 
ization of the interaction g rather than t. Our maximum 
estimate of S U]R < 0.23 meV would typically yield 5g/g < 
- 0.4 %, in agreement with the expected impact of the in- 
teractions termS. More specifically, if the interaction is 



assigned to the AF superexchange, we obtained an iso- 
topic change of Jaf compatible with the one, 5 J/ J ~ - 
0.6 %£, deduced from the Neel temperature in undoped 
cuprates. Finally, this simple analysis of the tiny shift (if 
any) of the resonance peak energy in term of the spin- 
exciton model shows that the various microscopic terms 
entering in its expression do not exhibit a strong isotope 
effect, thus placing a severe limit on the role of electron- 
phonon coupling in high-T c cuprates. 

In conclusion, using inelastic neutron scattering experi- 
ments, we observe no significant shift of the magnetic res- 
onance peak energy in YBCO upon substitution of oxy- 
gen 16 O by its isotope 18 O. In contrast to previous claims, 
this suggests that the spin-phonon coupling is negligible 
in high-T c cuprates near optimal doping. The absence of 
a measurable effect on the INS resonance mode does not 
however exclude the possibility that isotope substitution 
can have a noticeable effect on the magnetic properties 
at much lower doping. 

We thank Ilya Eremin and M.V. Eremin for stimulat- 
ing discussions, alerting us to possible isotope effects in 
the resonance energy. 
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